Relaxation dynamics in strained fiber bundles 
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Under an applied external load the global load-sharing fiber bundle model, with individual fiber 
strength thresholds sampled randomly from a probability distribution, will relax to an equilibrium 
state, or to complete bundle breakdown. The relaxation can be viewed as taking place in a 
sequence of steps. In the first step all fibers weaker than the applied stress fail. As the total load 
is redistributed on the surviving fibers, a group of secondary fiber failures occur, etc. For a bundle 
with a finite number of fibers the process stops after a finite number of steps, t. By simulation and 
theoretical estimates, it is determined how t depends upon the stress, the initial load per fiber, 
both for subcritical and supercritical stress. The two-sided critical divergence is characterized by 
an exponent —1/2, independent of the probability distribution of the fiber thresholds. 



PACS numbers: 62.20.Mk 



INTRODUCTION 



the applied external load per fiber, 



Bundles of fibers, with a statistical distribution of 
breakdown thresholds for the individual fibers, are simple 
and interesting models of failure processes in materials. 
They can be analyzed to an extent that is not possible 
for most materials (for reviews, see [H, 0, H, HI IH). 

We consider a bundle with a large number N of elas- 
tic and parallel fibers, clamped at both ends. When the 
load on fiber i is increased beyond a threshold value x», 
the fiber ruptures. The breakdown thresholds Xj for the 
separate fibers are assumed to be independent random 
variables with a probability density p(x), and a corre- 
sponding cumulative distribution function P(x): 



Prob (xi < x) = P(x) 



P{y) dy. 



(1) 



The mechanism for how the extra stress caused by a fiber 
failure is redistributed among the unbroken fibers must 
be specified. We study here the classical version, the 
equal-load-sharing model, in which a ruptured fiber car- 
ries no load, and the increased stress caused by a failed el- 
ement is shared equally by all the remaining intact fibers 
in the bundle [f}]. 

If an external load F is applied to a fiber bundle, the 
resulting failure events can be seen as a sequential process 
0, H, [§] ■ In the first step all fibers that cannot withstand 
the applied load break. Then the stress is redistributed 
on the surviving fibers, which compels further fibers to 
fail, etc. This iterative process continues until all fibers 
fail, or an equilibrium situation with a nonzero bundle 
strength is reached. Since the number of fibers is finite, 
the number of steps, t, in this sequential process is finite. 
In this paper we determine how t depends upon the num- 
ber of fibers and, more importantly, upon the stress ct, 
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a = F/N. 



(2) 



At a force x per surviving fiber, the total force on the 
bundle is x times the number of intact fibers. The ex- 
pected or average force at this stage is therefore 



F(x) =Nx(l- P(x)). 



(3) 



One may consider x to represent the elongation of the 
bundle, with the elasticity constant set equal to unity. 
The maximum F c of F(x) corresponds to the value x c 
for which dF/dx vanishes. Thus 



1 - P(x c ) - x c p(x c ) = 0. 



(4) 



We characterize the state of the bundle as subcritical or 
supercritical depending upon the stress value relative to 
the critical stress 



Fc/N, 



(5) 



above which the bundle collapses completely. Critical 
properties of fiber bundles have been discussed before, 
but with a signature different from the one that we use 
here, and always with the critical point approached from 
the subcritical side 0, H, [13]. The function t(a) that 
we focus on, however, exhibits critical divergence when 
the critical point is approached from either side. As an 
example, we show in Fig. 1 t(a) obtained by simulation 
for a uniform threshold distribution. 
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FIG. 1. Number of relaxation steps t(cr) for a fiber bundle with 
a uniform threshold distribution {TTJ. Here <r c = 0.25. The figure 
is based on 1000 samples, each with N = 10 6 fibers. 

We study the stepwise failure process in the bundle, 
when a fixed external load F — Na is applied. Let N t be 
the number of intact fibers at step no. t, with Nq = N. 
We want to determine how Nt decreases until the degra- 
dation process stops. With N t intact fibers, an expected 
number 



[NP(Na/N t )} 



(6) 



of fibers will have thresholds that cannot withstand the 
load, and consequently these fibers break immediately. 
Here [X] denotes the largest integer not exceeding X. 
The number of intact fibers in the next step is therefore 



N t+1 =N- [NP(N*/N t )} 
Since N is a large number, the ratio 



Nt 
N 



(8) 



can for most purposes be considered a continuous vari- 
able. By © we have essentially @, H, [§] 



nt+i = 1 - P{a/n t ). 



(9) 



In Sec. II we study t{a) in the supercritical domain, 
while Sec. Ill is devoted to subcritical situations. Simula- 
tion results are presented for two threshold distributions, 
the uniform and a Weibull distribution, and these are 
compared with detailed analytic results. The theoretical 
analysis is, however, not limited to these special thresh- 
old distributions. In Sec. IV we summarize our results 
and discuss briefly the approximations involved. 



II. SUPERCRITICAL RELAXATION 

We investigate first the supercritical situation, a > a c , 
with positive values of 



e = a-a n 



(10) 



and start with the simplest model. 



A. Uniform threshold distribution 

Consider a bundle in which the failure thresholds are 
distributed according to the uniform distribution 



P(x) 



x for < x < 1 
for x > 1 



(11) 



For this case the load curve is parabolic, 

F = Nx(l — x), (12) 

with the critical point at x c — 1/2, a c — 1/4. Simu- 
lation results for the uniform threshold distribution are 
presented in Fig. 1. 

The basic equation J9]) takes the form 



nt+i 



1-^ 

n t 



1 



rii 



(13) 



This nonlinear iteration can be transformed into a linear 
one by the following procedure. Introduce first 



into (fl3"ll . with a result that may be written 
Vt+i ~ Vt 



i + ytyt+i 



(7) Putting 



y t = t&nvt, 



(14) 



(15) 



(16) 



we have 
2^1 



tani> t _|_i — tanz; t 



1 + tani'i+i tan t>4 
Hence vt+i — v t = tan _1 (2y / e), with solution 

v t = v + t tan _1 (2 v / e). 
In the original variable the solution reads 



tan(v t+ i - v t ). (17) 



(18) 



n t 



f tan tan 



i 

-1/2 



no , 



+ t tan- 1 (2 v / i) )(19) 



= i-y^tan(-tan 1 (l/2s/e) + t tan 1 (2 v / ^))(20) 

where n = 1 has been used. 

Eq. (fl3|) shows that when n t obtains a value in the 
interval (0,<r), the next iteration gives complete bundle 
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failure. Taking n t = a as the penultimate value gives a 
lower bound, ti, for the number of iterations, while using 
nt = in ((20|) gives an upper bound t u . Adding unity 
for the final iteration, (|20|) gives the bounds 



tu(o-) = 1 + 



2 tan" 1 (1/2 Ve) 
tan- 1 (2Ve) 



and 



1 



tan-^Cl - f)A/f) + tan-^l^yg) 
tan" 1 (2Ve) 



(21) 



(22) 



It is easy to show that t u (<r) — U{a) — 1. In Fig. 2 A we 
show that these bounds nicely embrace the simulation 
results. 
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FIG. 2. Simulation results with supercritical stress for (A) 
the uniform threshold distribution JTTJ, and (B) the Weibull 
distribution J24j. The graphs are based on 10000 samples with 
N = 10 6 fibers in each bundle. The dashed lines represent the 
theoretical estimates J2TJ, lt22l and J32J. 

Note that both the upper and the lower bound behave 
as e~2 for small e. A rough approximation near the crit- 
ical point is 



t{a) w ±7r(er - a c ) I. 



(23) 



B. General threshold distributions 

The uniform distribution is amenable to analysis to 
a degree not shared by other threshold distributions. 
Therefore we now discuss how to handle other distribu- 
tions, and we start by a special Weibull distribu- 
tion of index 5, 



P(x) = l-e 



x>0. 



(24) 



The critical parameters for this case are x c — 5 1 / 5 = 
0.72478 and a c = {be)- 1 ' 5 = 0.5933994. Simulation re- 
sults for t(a) are displayed in Fig. 2B for the Weibull 
supercritical case. The variation with the external stress 
a is qualitatively similar to the results for the uniform 
threshold distribution. 

The interesting values of the external stress are close 
to <7 C , because for large supercritical stresses the bundle 
breaks down almost immediately. For a slightly above a c 
the iteration function 



n t +i = /(nt) = 1 - P(<t/th) = e^/ n ^ , 
takes the form sketched in Fig. 3. 



(25) 




FIG. 3. The iteration function f(n) for the Weibull distribution 
||24|I . Here a = 0.6, slightly greater than the critical value 
cr c = 0.5933994. 

The iteration function is almost tangent to the reflec- 
tion line n-t+i = n t and a long channel of width propor- 
tional to e appears. The dominating number of iterations 
occur within this channel (see Fig. 3). The channel wall 
formed by the iteration function is almost parabolic and 
is well approximated by a second-order expression 

n t +i = n c + {n t - n c ) + a(n t - n c ) 2 + b(a c - a). (26) 

Here n c = e -1 / 5 is the fixed point, n t +i — nt, of the 
iteration at a = a c . With u — (n — n c )/b and e = a — a c 
(l26l) takes the form 



u t+ i - u t 



-Auf 



(27) 
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with A = ab. In the channel u changes very slowly, so 
we may treat the difference equation as a a differential 
equation: 



du , 9 
— = -Au 2 
dt 



with solution 

t\f~Ae — — tan" 



constant . 



(28) 



(29) 



Thus 



tf-U = {Ae)-i { 



tan 



tan" 



(30) 

is the number of iterations in the channel, starting with 
Ui, ending with Uf. This treatment is general and can 
be applied to any threshold distribution near criticality. 
Although the vast majority of the iterations occur in the 
channel, there are a few iterations at the entrance and at 
the exit of the channel that may require attention in spe- 
cial cases. The situation is similar to type I intermittency 
in dynamical systems [ll|, but in our case the channel is 
traversed merely once. 

For the Weibull distribution the expansion (|26| has the 
precise form 



-5 1/5 
2 e 



" ~ e~ 1/5 + (n - n c ) 

(n-n c ) 2 -^^{a-ac), (31) 



where n c 



|e 1 /5 j 5 = 5 i/5 andA=|(5e) 



1/5 



For completeness we must also consider the number 
of iteration to reach the entrance to the channel. It 
is not meaningful to use the quadratic approximation 
([31"]) where it is not monotonously increasing, i.e. for 



n > n m = n c + l/(2a) = fe" 1 / 5 ~ 0.98. Thus we 
take n, = n m as the entrance to the channel, and add 
one extra iteration to arrive from no = 1 to the channel 
entrance. (Numerical evidence for this extra step: For 
a = a c the iteration (|25l) starts as follows: no = 1.00, 
ni = 0.93, «2 = 0.90, while using the quadratic function 
with n — n m = 0.98 as the initial value, we get after 
one step m = 0.90, approximately the same value that 
the exact iteration reaches after two steps.) With nf = 
we obtain from lf30l in the Weibull case the estimate 



t = l + (Ae)- 1/2 {tan^ 1 (e- 1/5 v /I7^/56) 

+ tan- 1 (e~ 1/5 VA 7 e/fr)}, (32) 

with A = f (5e) 1/5 and b = 5 1 / 5 . 

We see in Fig. 2B that the theoretical estimate i|32|) 
gives an excellent representation of the simulation data. 
Near the critical point l132[) has the asymptotic form 



t m ir(Ae)- 1/2 oc (a - a c y 1/2 , 



(33) 



with the same critical index as for the uniform threshold 
distribution. The divergence is caused by the large num- 
ber of iterations in the narrow channel in Fig. 3. For a 



general threshold distribution such a channel will always 
be present, and therefore the divergence (|33|) is univer- 
sal. Moreover, the amplitude of (er — a c Y 1 ' 2 , as well as 
an excellent representation of the complete t{a) function 
may, for a given threshold distribution, be obtained by 
the same procedure as used above for the Weibull case. 



III. SUBCRITICAL RELAXATION 

We now assume the external stress to be subcritical, 
a < a c , and introduce the positive parameter 



e = Or - er 



(34) 



to characterize the deviation from the critical point. 

Also in the subcritical situation a bundle with the uni- 
form threshold distribution ifTTj) can be analyzed analyti- 
cally to a greater extent than for other distributions, and 
consequently we start with this case. 



A. Uniform threshold distribution 

Using a similar method as in the supercritical situation 
we introduce into 11131) 



n t = l + \fejz u 
as well as o = \ — e, with the result 



2^- 



Zt+i - z t 
1 - z t +i z t ' 

In this case 

z t = tanh Wt 
is the useful substitution. It gives 



(35) 



(36) 



(37) 



tanh Wt+i - tanh w t 

= — : — - — = tanh(w t+ i - w t ). (38) 

1 — tannu! t+ i tanniy t 

Thus wt+i — w t — tanh _1 (2v / e), i.e. 

w t = w + t tanh _1 (2Ve). (39) 
Starting with no = 1, we obtain zq = 2-y/e and hence 

Wt = (1 + t) tanh _1 (2Ve). (40) 
This corresponds to 



n t 



tanh {(1 + t) tanh -1 (2^)} 



(41) 



in the original variable. 

Apparently n t reaches a fixed point n* = i + ^fe after 
an infinite number of iterations. However, our bundle 
contains a finite number of fibers, and therefore only a 



5 



finite number of steps is needed for the iteration to arrive 
at a fixed point N* of the integer iteration J7]) , 



N t+1 =N- [aN 2 /N t ] 



(42) 



Since X < [X] < X + 1 a fixed point N* of (@2l) must 
satisfy 

y (1 + VI - 4cr) < AT* < ~ (V + 1 + v /A r2 (l-4cr)+2Af + l 

(43) 

It is interesting to note that l]42l) has in general several 
fixed points for a given value of a. With N — 10 6 and 
a = 0.249, for instance, there are nine fixed points, viz. 
531623, 531624, . . . , 531631, the complete set of integers 
within the interval (|43|) . Since our iteration starts high 
at iV = N, with steadily decreasing values of N t , it 
will stop at the upper fixed point, the largest integer 
satisfying (|43"1) . 

As long as N(l — 4<j) ^> 1, which is fulfilled in our 
simulations, we may take 



N* 



N 



l-VT^4^) + i(l + (l-4a)- 1 / 2 ) (44) 



as a good approximation to the upper fixed point (in 
the example above lf44j) gives 531631.1, compared with 
N* = 531631). 

As a consequence we use 



N 2 v AN 1 



(45) 
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as the final value in lj41j) . Consequently we obtain the 
following estimate for the number of iterations to reach 
this value: 

t(a) = -l + C ° th ~ 1{1 + ^ +2 ^ 4N£ \ (46) 
tanh" :L (2Vi) V ' 

We see in Fig. 4A that the simulation data are well ap- 
proximated by the analytic formula (|46j) . 
For very large N |46|) is approximated by 



MgjVg) -i/2 

4 



(47) 



The critical behavior is again characterized by a square 
root divergence, in this case somewhat modified by a log- 
arithmic term. 



B. General threshold distribution 



Again we use the Weibull distribution l|24j) as an ex- 
ample threshold distribution. Simulation results for the 
subcritical Weibull distribution are shown in Fig. 4B. 



FIG. 4. Simulation results with subcritical stress for (A) the uni- 
form threshold distribution, and (B) the Weibull distribution l|24|l . 
The graphs are based on 10000 samples with N = 10 6 fibers in each 
bundle. The dotted lines are the theoretical estimates, eq. II46II for 
the uniform distribution case, and eqs. I|55I56|I for the Weibull case. 

Forgetting for the moment the finiteness of the fiber 
bundle, the iteration |(9]), 



nt+i = 1 - P{a/n t ), 



(48) 



will reach a fixed point n* after infinite many steps. The 
deviation from the fixed point, nt — n* will decrease ex- 
ponentially near the fixed point 0, @, S E2] : 



nt — n oc e 



-t/r 



(49) 



with 



t = l/]n{n* 2 a- x /p(a/n*)} . (50) 
For the Weibull threshold distribution, in particular, 

p(a/n*) = 5(cr/n*) 4 exp (->/n*) 5 ) = 5a 4 /n* 3 , (51) 
and thus 

r = l/ln(n* 5 /5cr 5 ) (52) 
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for the Weibull case. If we allow ourselves to use the 
exponential formula lj49j) all the way from no = 1, we 
obtain 



(1 -n*)e-'/ T . 



(53) 



For a finite number N of fibers the iteration will stop 
after a finite number of steps. It is a reasonable suppo- 
sition to assume that the iteration stops when N t — N* 
is of the order 1. This corresponds to take the left-hand 
side of lj53j) equal to 1/N, The corresponding number of 
iterations is then given by 



t = t ln(JV(l - n*)) 



In (N(l - n*)) 
ln(n* 5 /5cr 5 ) 



(54) 



Solving the Weibull iteration n* = exp(— (a/n*) 5 ) with 
respect to a and inserting into lf54"l) . we obtain 



* = 



ln{iV(l-n*)} 
~ ln{5(-lnn*)} 
n^-lnn*) 1 / 5 . 



(55) 
(56) 



These two equations represent the function t{a) on pa- 
rameter form, with n* running from n c = e -1 ^ 5 to 
n* = 1. In Fig. 4B this theoretical estimate is compared 
with the simulation data. The agreement is satisfactory. 



For n* 



-1/5 



55j) shows that t is infinite, as 



it should be. To investigate the critical neighborhood we 
put n* = rt c (l + £;) with £ small, to obtain to lowest order 



t = oic 1 ) 

This gives, once more, the divergence 
t(a) oc (a c - a)- 1 ' 2 . 



(57) 
(58) 



(59) 



For a general threshold distribution the same proce- 
dure may be followed. To use the exponential approach 



to the fixed point, as we have done, may seem to be 
doubtful. But the rational is that for small a the start- 
ing point n — 1 is already rather close to the fixed point, 
while for larger a it doesn't matter much if the first few 
iterations is not described well by the exponential for- 
mula, since in any case the majority of the iterations 
occur close to the fixed point. 



IV. CONCLUDING REMARKS 

A detailed numerical and analytic study of the relax- 
ation dynamics in finite fiber bundles subjected to ex- 
ternal loads is presented. The relaxation takes place in 
a number, t(a), of steps: In each step all fibers weaker 
than the load per surviving fiber burst, and the relax- 
ation proceeds until equilibrium is reached or until all 
fibers have failed. As function of the initial stress a the 

1/2 

number of steps, t(a), shows a divergence |cr — cr c | ' 
at the critical point, both on the subcritical and super- 
critical side. This is a generic result, valid for a general 
probability distribution of the individual fiber breakdown 
thresholds. 

On the supercritical side t(a) for large N, is indepen- 
dent of the system size N. On the subcritical side there 
is, however, a weak (logarithmic) A-dependence, as wit- 
nessed by eqs. (46), (47) and (55). 

The analytic estimates are based on the average 
strength of a group of fibers. The comparison with the 
simulation data on the probabilistic distribution of fiber 
strength shows that this is a satisfactory calculation pro- 
cedure. 
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